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ABSTRACT 

Feynman rules for the vacuum amplitude of fermions coupled to external gauge and 
Higgs fields in a domain wall lattice model are derived using time-dependent perturbation 
theory. They have a clear and simple structure corresponding to 1-loop vacuum graphs. 
Their continuum approximations are extracted by isolating the infrared singularities and 
it is shown that, in each order, they reduce to vacuum contributions for chiral fermions. In 
this sense the lattice model is seen to constitute a valid regularization of the continuum 
theory of chiral fermions coupled to weak and slowly varying gauge and Higgs fields. 
The overlap amplitude, while not gauge invariant, exhibits a well defined (modulo phase 
conventions) response to gauge transformations of the background fields. This response 
reduces in the continuum limit to the expected chiral anomaly, independently of the phase 
conventions. 
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1 Introduction 



Recent efforts to provide a lattice regularized model for chiral fermions have made some 
progress using Kaplan's idea of the domain wall [1]. This can be formulated as a 4- 
dimensional Euclidean lattice embedded in 4+1-dimensional spacetime. The continuous 
time coordinate is an artificial variable whose purpose is to accomodate a barrier, or 
domain wall, with which 4-dimensional chiral fermions can be associated. The barrier 
corresponds to a mass term discontinuity in the time-dependent Hamiltonian of the 4+1- 
dimensional system. Couplings to time-independent (external) gauge and Higgs fields are 
included in this Hamiltonian. The aim is to compute the vacuum transition amplitude 
for this 4+1-dimensional system and extract from it a quantity that can be interpreted 
as the Euclidean vacuum amplitude for chiral fermions on a 4-dimensional lattice. 

An equivalent picture, developed by Narayanan and Neuberger [2,3], interprets the 
4-dimensional Euclidean amplitude as the overlap of ground states belonging to the two 
distinct Hamiltonians that govern the 4+1-dimensional system on either side of the bar- 
rier. 

The motivation behind these efforts is to obtain a model suitable for numerical, i.e., 
non-perturbative, studies of chiral theory. It is not yet clear that this aim will be achieved, 
but some encouraging results have been obtained by Narayanan and Neuberger [3,4] in 
the context of a 2-dimensional model whose continuum version is soluble. 

Our purpose here is more limited. We want to show that the domain wall or overlap 
prescription is perturbatively correct. By this we mean that, for weak and slowly varying 
external gauge and Higgs fields, the prescription yields, in each order, a regulated version 
of the Feynman graph that one expects to find in a continuum theory of chiral fermions. 
Some work in this direction is already available in the literature. On the one hand, a 
set of rules for computing the perturbative contributions to the overlap amplitude was 
developed by Neuberger and Narayanan [3]. On the other hand, the continuum limit of 
low order contributions has been examined by several groups [2,5,6]. The latter work is 
concerned mainly with establishing that the expected chiral anomalies are indeed present. 
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Recently, however, we verified that the vacuum polarization tensor for chiral fermions of 
the SU(2) x U{\) standard model is recovered in the continuum limit [7]. Kaplan and 
Schmaltz [8] have shown in the continuum version, that the phase of the overlap coincides 
with the 7/-invariant of Alvarez-Gaume et al. [9]. This is a non-perturbative result. Some 
non-perturbative analytic work in 2-dimensions is also available [10]. 

To set up the lattice model that corresponds to a collection of Weyl fermions coupled to 
background gauge and Higgs fields in a 4-dimensional Euclidean spacetime, one begins by 
doubling the number of fermion components. Corresponding to each 2-component Weyl 
fermion, or ip R (x), introduce a 4-component Dirac field, ip(n,t), defined on the 

sites of a 4-dimensional integer lattice, G TL A . The "time" coordinate, t, is continuous. 
Construct the Hamiltonian as a bilinear form, 

H{t) = ]T ip(n, ty H(n, m, t) ip(m, t) (1.1) 

n,m 

where the matrices H(n,m,t) are covariant functionals of the background fields. These 
background fields, A^(x) and (f)(x), are assumed to be smooth and independent of t. The 
functional dependence of the matrices H(n, m, t) is specified in detail in Appendix A. 
Here, we remark only that their time-dependence is confined to the mass-like term, 

e(t)A l5 T c 5 nm (1.2) 

where e(t) = sign(t), A is a positive parameter representing the height of the domain wall 
and T c is a diagonal matrix with eigenvalues +1(— 1) corresponding to right(left) handed 
Weyl fermions. It commutes with the Dirac algebra. See Appendix A for details. 

The structure of the Hamiltonian (1.1) is very simple. It has a discontinuity at t — 
but is otherwise independent of time, 

f H + (A), t > 
H(t) = (1.3) 
[ H_(A), t<0 

where the argument, A, stands collectively for the background gauge and Higgs fields. 
Although the Hamiltonian is discontinuous, the Heisenberg-picture field, ip(n,t), is con- 
tinuous at t — where it coincides with the Schroedinger-picture field, ip{n). When the 
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background fields are weak, as we shall assume, there is a natural separation into free and 
interaction terms, 



and one can set up the usual perturbation series. The only unusual feature here is the 
discontinuity in the free Hamiltonian. 

With the two Hamiltonians, H + (A) and H_(A), one can construct two distinct normal- 
ized ground states, \A+) and \A—), respectively. Of particular interest is the functional 
T(A) defined by the overlap, 



where |±) denotes the respective ground states of the free Hamiltonians, H ± = H±(0). 
This is the functional whose perturbation development we shall consider and which, we 
shall show, reduces in the continuum limit to the connected vacuum amplitude for a set 
of Euclidean Weyl fermions. An efficient way to compute T(A) is by means of time- 
dependent perturbation theory. The rules for expressing the contributions in terms of 
connected 1-loop vacuum graphs are obtained in Sec. 2. 

Since the free Hamiltonian depends on time, the free particle Green's function will 
not be invariant under time translations and, as a result, the detailed structure of the 
perturbative contributions to T(A) is more complicated than in familiar theories. However, 
these complications tend to become unimportant in the continuum limit. When the 4- 
momenta carried by external fields are small compared to barrier height and inverse lattice 
spacing, amplitudes are dominated by infrared singularities, i.e. thresholds associated 
with the propagation of light fermions. The leading infrared singularities are insensitive 
to lattice structure and can be computed by an effective continuum field theory of chiral 
fermions in 4-dimensional Euclidean spacetime. It is precisely this continuum theory that 
the lattice model regulates. The emergence of a chiral continuum theory is discussed in 
Sec. 4. It depends crucially on the infrared behaviour of the free fermion Green's function 
whose detailed structure is considered in Appendix B. 



H(t) = H (t) + V 



(1.4) 




(1.5) 
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The functional T(A) is not gauge invariant. Under a gauge transformation of the 
background fields, A — > A e , it responds according to 

T(A e ) = r(A) - i $+(0, A) + i A) (1.6) 

where $± are real angles associated with transformations of the ground states, \A±). 
There is some arbitrariness in these angles that reflects the role of phase conventions 
in the construction of \A±). One possibility would be to impose the Brillouin-Wigner 
convention: that the overlaps (±|A±) shall be real and positive. This choice was adopted 
by Narayanan and Neuberger in their original formulation of the overlap prescription 
[2]. It is well adapted to time-independent perturbation theory and was used also in 
our low order computations [5,7,11]. Here we shall use another convention that is better 
adapted to time-dependent perturbation theory, to be explained in Sec. 2. With either of 
these conventions the difference, $ + — is non-vanishing in general. In both of them, 
however, it can be shown that $ + — $_ reduces to the standard chiral anomaly in the 
continuum limit. This means, in particular, that T(A) becomes gauge invariant in the 
continuum limit if the Weyl fermions belong to an anomaly-free combination [5]. 

Gauge transformations are discussed in Sec. 3 and the angles <3>± are defined there. The 
relative phase between the ground states of this paper and the Brillouin-Wigner states 
used in earlier work is discussed in Appendix C where the time-independent formalism is 
briefly reviewed. 

A subtle point concerning gauge transformations and the continuum limit is raised 
in Sec. 4. This limit exhibits a lack of "uniformity". One finds that the gauge variation 
of the continuum limit of the effective action differs from the continuum limit of the 
gauge variation, $ + — This is because the continuum limit of T(A) is dominated 
by infrared singularities that are not present in $ + — The latter quantity, it will be 
seen, is determined by massive fermions and reduces to a local form, the integral over 
4-dimensional Euclidean spacetime of a pseudoscalar density in the slowly varying gauge 
fields and their derivatives. The coefficients in this density are finite lattice-dependent 
quantities. The chiral anomaly is a sub-dominant effect from the infrared point of view. 
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(This is only to be expected since, in continuum gauge theory the anomaly arises in a 
parity violating amplitude which is ultraviolet convergent and unambiguous, but whose 
gauge variation is, at least superficially, ultraviolet divergent.) In order to recover the 
standard consistent anomaly of continuum gauge theory from the local expression for 
$ + — $_ it is necessary to let the barrier height, A, become vanishingly small relative to 
the lattice cutoff, a -1 . The condition, Aa <C 1, was used in Ref. [5] where we obtained 
the chiral anomaly by computing $ + — This condition is used implicitly in many 
studies of the overlap prescription - so called continuum models - where the fermions 
are represented by smooth fields in 5-dimensional spacetime, and ultraviolet questions 
are ignored [2,3,11,12,8]. However, it should be recognized as a non-essential technicality. 
In Sec. 4 and Appendix B we show that the continuum theory emerges as the infrared 
dominant part of T(A) provided only that the background fields are slowly varying on the 
scale of A^ 1 . It is not necessary to assume A -1 ^> a.[] 

In this paper we are exclusively concerned with the domain-wall-overlap formulation 
of chiral gauge theories on the lattice. For some other approaches see [13-17]. 

2 Perturbation theory 

Since the domain wall problem is unusual in having a time dependent free Hamiltonian 
we begin with a brief description of perturbation theory in the interaction picture. We 
assume that the free Hamiltonian is invariant with respect to lattice translations so that 
Fourier transforms can be used. The interaction picture equations of motion take the 
form 



: We wish to acknowledge a useful correspondence with H. Neuberger who insisted on this point. 



i dt ip(p,t) 



[ip{p,t),H (t)] 



H(p,t)i>(p,t) 



(2.1) 
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where H(p,t) is an hermitian matrix, discontinuous at t — but otherwise independent 
of time (see Appendix A). The solution of (2.1), continuous at t — 0, is given by 

f e -itH+(p) t>Q 

Hp,t) = \ (2.2) 
[ e~ itH -^ rj;(p), t<0 

The time-dependent states of the interaction picture are governed by the usual unitary 

operator, 

'Tfe-'fo^'A, t>0 
n(t) = { ) { (2.3) 

where T denotes antichronological ordering. In order that these integrals converge for 
t — > ±oo, the operator V(t) should include the damping factor, e" e '*', i.e. 

% d t V(t) = [V(t), H (t)] - %e sgnit) V(t) (2.4) 

There are two free fermion ground states, |+) and |— ), defined by 

H 0+ \+} = 0, #o-|->= (2.5) 

where the respective Dirac seas are filled. Expressions for the 1-body Hamiltonians H±(p) 
are given in Appendix A. The details are not important for now, except for the existence 
of a gap, 

\H±(p)\ > k^min 

(A)>0 

where A represents the height of the domain wall. Ground states of the interacting theory 
can be generated adiabatically from these states [18]. Thus, for small but finite e define 
the asymptotic states 

\in ±) = IX-oor^i) 
\out±) = IXoor^i) ( 2 - 6 ) 

When e tends to zero these states converge, apart from a singular phase, onto eigenstates 
of the Schroedinger-picture Hamiltonians, 

H ± (A) = H 0± + V (2.7) 



where A denotes a collection of time-independent external fields. One can show [18] 

lim (\in ±) e"^ ±/l£ ) = \A±) 

lim (\out ±) e^ ±/i£ ) = \A±) (2.8) 

where the functional E±(A) are related to the ground state energies defined by 

H ± (A)\A±) = \A±) E ± {A) (2.9) 

The relation between E and E takes a simple form when both are expanded in powers of 
V, viz. 

E = Ei + E 2 + E 3 + . . . 

E = E ± + ^E 2 + ^E 3 + ... (2.10) 

To compute these energies using time-dependent perturbation theory one writes, for 
example, 

(out + \in+) = e 2 ^ +/i£ 

= (+\n(oo) ^(-oo)- 1 ^) 

= (+\T(e-*^ dtVit) )\+) 
which can be reduced to the computation of connected vacuum graphs, 

E + = lim | (+|T (e-/^ 1 v ®) \+) con (2.11) 
Since the interaction is bilinear, 

V(t) = J {^j (^j V(p,*) f V(p,q) Mq,t) e- £ l*l (2.12) 
There is only one connected graph in each order. One finds, 

E + = ^4f- l°° dti ...dt N e- e(l * l|+ - +|tjv|) • 

2 jy iV J — CO 

7 (l 1 ) 4 '--^) 4 tr[v(p uP2 )S + (p 2 M-h)V(p 2 ,p,)... 

...V(p N , Pl )S + ( Pl ,t N -t 1 )} (2.13) 



where the propagator S+ is defined by 

(+m(q,t) V(p,O f )l+) = (27r) 4 M(z-p) S+(q,t-1?) (2.14) 

It is invariant with respect to time translations because the time dependence of ip(q,t) is 
determined here by the time-independent Hamiltonian, H 0+ . Indeed, we can write 

S + (q, t-t')= r — - TT , / — - e~ iE ^ (2.15) 

+ w ' ; J-oo 2?r E-H+(q) + irj sgn(H + ) v ; 

with i] > 0. The momentum integrations in (2.13) are compact. They range over the 
Brillouin zone, a cell of volume (27r) 4 in lattice units. The delta function in (2.14) is 
periodic with respect to reciprocal lattice translations. 

An expression analogous to (2.13) can be written for E_ by replacing S + with S-. 
More interesting is the transition amplitude that includes the barrier effect, 

(out+\in-) = (+\T (e-'f-o** V{t) ^j |-> (2.16) 

where the time dependence of tp(q,t) is now determined by the time-dependent Hamilto- 
nian, H (t). Here one computes the connected vacuum graphs using the propagator Sf, 
defined by 

(+\T(^(q,t) V(p,O f )|-> = (+h) (2vr) 4 Mg - p) S F (q,t,t') (2.17) 

This propagator is not invariant with respect to time translations. It is very much more 
complicated than (2.15) and it includes long range effects due to the propagation of chiral 
fermions. Its detailed structure is discussed in Appendix B. 

The order iV contribution to —(out + \in—) con is given by an expression analogous to 
(2.13), 

/ .\ AT 

fdt 1 ...dt N e -e(N+-+l*"l) • 



:-t) N 



N 



/(f) ■■■( 



dp N \ 
2tt / 



V(pi,P2)S F (p 2 ,t 1 ,t 2 )V(p 2 ,P3) . . . 

...V(p N ,p 1 )S F ( Pl ,t N ,t 1 )} (2.18) 
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Finally, these results can be put together to give the effective action functional, T(A), 
defined by (1.5), 

(+|_) e -r(A) = (A + \A-) 

= (out + \in-) e -^+ + *-V i£ 

= (out + \in—)(out+ \in+)^ 1 ^ 2 (out — \in— }~ 1 / 2 

or, in terms of connected vacuum graphs, 

F(A) = -(out + \in-) con + ]- (out + \in+) con + ^ (out - \in-) con (2.19) 

The respective terms are to be computed using the propagators, Sp,S + and 5_. The 
limit e — > is understood. In fact, the auxiliary pieces, (out ± \in±) con , are needed only 
to cancel the singularity in (out + \in— ) con as e tends to zero. The effective action is 
obtained as the regular part of —(out + \in—) con . 

3 Gauge transformations 

Infinitesimal time-independent gauge transformations are generated by the operator 

F e = S V(n) f 9(n) ip(n) 

= /(^) 4 (^) 4 m f o( P - q )m (3.i) 

where 8(n) is an hermitian matrix belonging to the algebra of the gauge group. It is 
slowly varying on the lattice and will be interpolated by a smooth function, 9(x), that 
defines the transformations of the background fields, A — > A 6 , such that 

e iFe H±(A) e- tFe = H±(A e ) (3.2) 

Since the ground states are non-degenerate, at least in perturbation theory, it follows 
that they must transform according to 

e tFe \A±) = \A e ±) e l ^ e ' A) (3.3) 
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where the angles $± are real. These angles provide a representation of the group. Thus, 
if the product of two group elements is defined by 

then it is easy to show that the corresponding composition rule for $± is given by 

$(0 1 ,A e *) + $(0 2 ,A) = $(0 12 ,A) (3.4) 

Gauge transformations of the effective action, T(A), are obtained by substituting (3.3) 
into the definition (1.5), 

T(A e ) = F(A) - i $+(0, A) + % $_(0, A) (3.5) 

There is no reason to expect the difference, $ + — to vanish in general but we should 
expect some simplifications to occur when the background fields are slowly varying. To 
see what happens it is necessary to compute these angles in perturbation theory. 
To first order in 9 the transformation rule (3.3) takes the form 

8o\A±)=i{F -Q ± )\A±) 

which implies, for example, 

+iS e en(+\A+) (3.6) 

and likewise for This can be calculated in time-dependent perturbation theory using 
the adiabatic formula (2.8). Firstly, 

i 5 9 £n(+\A+) = i 5 e £n(+\in+) - ^ 5 e E + 



f o (+\T(5eV(t)e-^ dt ' v ^)\+) , 

/ dt " / — T v " 

J-oo (+\in+) e 



-oo (+\in+) 
A suitable formula for SeV(t) can be extracted from (3.2), 

S e V(t) = e itHo + SgV e- ttHo + e £t 

= e itHo+ i[F e , H 0+ + V] e~ itHo+ e £t 
= e £t t[F 6 (t),H 0+ ]+t[F 6 (t),V(t)] 

= -e £t d t Fg(t)+i[Fg(t),V(t)} 
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- - 5 e E + (3.7) 



which implies 

<+|T(W) e- i ^ dt ' V ^\+} = 

= -e £t d t (+\T fat) e-** —"^ |+) 

-{e £t - 1)(+|T (i[F e (t),V(t)]e^£-" v( n) |+> (3.8) 



(MT(F e (t)e-*f-~ dt ' v «A\ + ) 

im / dt < 1 <h 



For Abelian symmetries the result is relatively simple. In such cases the energy E + (A) 
is invariant in each order so that SgE + = 0. Also, the commutator [i*0, V] vanishes and 
the second part of (3.8) is absent. It follows that 

$ + = (+l^+) _ lim/ dte st di 

(+\m+) e^oj-oo (+\m+) 
= lim e fjt e £t (+\T fat) e f°- dt ' v ^ \+) con ( 3 .9) 

For non-Abelian symmetries the second part of (3.8) cannot be ignored. It contributes 
a singular term that cancels 5gE/e from (3.7) and a regular term that must be retained. 
Hence we can write 

= reg f_Jt{+\T [{e e*F e {t) - (e £ * - 1) i[F e (t), V(t)}} e~ l f-™ dt ' v(t '^ |+> (3.10) 

meaning the regular part at e = 0. A comparison of this result with the correspond- 
ing angle obtained in the time-independent formalism using the Brillouin-Wigner phase 
convention is made in Appendix C where (3.10) is computed up to second order in V(0). 

4 Infrared behaviour 

Having obtained general expressions for the perturbative contributions to the vacuum 
amplitude, {out + \in—), we now consider how to approximate them when the external 
fields are slowly varying. 

In each order, the formula (2.18) corresponds to a 1-loop vacuum graph constructed 
from the vertices, V(p, q), and propagators, Sp(p,t,t'). The vertex is expressible as an 
expansion in powers of the gauge field, A^(x), the first two terms of which are given by 
(A. 12). In general it should include powers of the Higgs field, (f)(x), as well, but we are 
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simplifying the structure by choosing to be constant and incorporating it in the free 
propagator, Sf- This propagator, discussed in Appendix B, is given by (B.14). While its 
general structure is quite complicated, near p = = it simplifies to the form given by 
(B.15), which exhibits a pole corresponding to the propagation of light chiral fermions. 

The loop integration in (2.18) ranges over the Brillouin zone, a cell of volume (27r) 4 
in lattice units. There is, of course, no ultraviolet divergence. What interests us here is 
the possibility of infrared divergences which we expect to dominate the amplitude when 
the gauge fields are slowly varying and the Higgs field is small. To see that there is an 
infrared singularity it is sufficient to examine the integrand of (2.18) in the vicinity of the 
point, pi — p 2 — ■ ■ ■ — Pn — 4> — 0. In this region the vertex (A. 12) reduces to 

V(p,q)~-i^A fl (p-q) + ... (4.1) 

and the propagator is dominated by the pole term in (B.15), 

Sf(p, t, f ) - 75 A e-"W+M + ... (4.2) 

2 ip + <p • 1 

Integration over the time coordinates, h, . . . ,tw, is trivial in this approximation since 

2A 1 



/OO 
-00 



e + 2ik i 

in the limit, e — > 0. Hence the integrand of (2.18) reduces to 



Tit \ 1 + ^ Tc 1 Tit \ 1 + 

4{Pl - P2) o ■ 2 , 7 nn ■ ■ ■ 4{PN ~ Pi) 



N 



2 ifc + • T 2 + 

r ~. . 1 -+■ 

tr 



4(h) 7: -71 U x I 7 nn ■ ■ ■ 4( k N)- 



(4.3) 



2 1(^-^)4-0. T'" /v ' 2 ip' + <f)-T_ 
where p\ — p — k±,p2 = p — k\ — k 2 , . . . and the external momenta are constrained to 
satisfy k x + k 2 + . . . + k N = 0. For > 4 there is clearly an infrared singularity at 
ki — . . . — k N — = 0, because the loop integration then diverges at p — 0. For 
N < 3, derivatives of order 3 — iV with respect to external momenta also diverge. These 
singularities are of course threshold effects associated with the propagation of light chiral 
fermions near their mass shell. The expression (4.3) is what one would expect to find in 
a 4-dimensional continuum theory described by the (Euclidean) Lagrangian density 

C = - i4 + • T) 1 + 75Tc V (4.4) 
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Our point is that (4.3) emerges as the dominant infrared effect in the lattice model. 
Subleading terms, down with respect to (4.3) by powers of ki, . . . , kjsr, <t> could be com- 
puted by improving the approximate formulae (4.1) and (4.2), but they would be lattice 
dependent and should be interpreted as scaling violations, irrelevant in the continuum 
approximation. In this sense the lattice overlap, or domain wall, prescription constitutes 
an ultraviolet regularization of the continuum system (4.4). 

To summarize our approach: we seek to isolate the contributions that are singular in 
the infrared. The leading singularity is lattice-independent, sub-leading and non-singular 
quantities are sensitive to the lattice and should not be computed. Lattice dependent 
quantities are either not relevant to the continuum theory or they can be incorporated in 
counterterms. A detailed exposition of this approach as applied to the vacuum polariza- 
tion tensor (JV = 2) is given in Ref.[7]. 

An important qualification should be made. The expression (4.3) represents the dom- 
inant infrared contribution only if the propagator Sp(p,t,t') has no other poles. The 
argument assumes that p = is the only point in the Brillouin zone where the propaga- 
tor is singular. It must be shown explicitly that there are no other such points, i.e. no 
doubling of fermions [19,20]. This matter is dealt with in Appendix B. 

It may be remarked that, since the infrared dominant and lattice independent contri- 
butions (4.3) coincide exactly with the continuum theory formulae, they must also carry 
the expected chiral anomalies. For example, with N = 3 and, for simplicity, = 0, the 
parity violating part of the amplitude is given by 



(4.5) 

where the loop integration is understood to comprise a small region around p — 0. The 
asymptotically dominant contribution to T 5 when fci, &2, k 3 tend to zero is obtainable, 
however, by extending the range of p to the entire IR 4 since the resulting integral is, in 
fact, ultraviolet convergent. (There is no 5*0(4) invariant, pseudoscalar local term of 
dimension 4 that could serve as a counterterm.) One may calculate this amplitude and 
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verify that it satisfies 

6e rs(A) = "24^ / d ' X ^ tr ( Tc 9k Ax ^ K 9 ) (4 - 6) 
to second order in A. 

ft is interesting to compare the resuft (4.6) with the general formula (3.5) or, to first 
order in 9, 

5 e r(A) = -i<S> + (0,A)+i<f>_(e,A) (4.7) 

The angle <3> + is given by (3.10) which can be expanded in powers of the interaction, 
V. The result up to terms of second order is given by (C.14). In the time-dependent 
formalism the angles <£>+ and <£>_ are computed using the Green's functions S + and S-, 
respectively. These propagators, discussed in Appendix B, do not have a long range 
structure. They are regular at p — 0, 

S ± (p,t-t')= 1 -(e(t-t')± l5 T c ) e-^-*'l + o(|) (4.8) 

This means that the functionals $ + and $_ do not have any singularities. For slowly 
varying fields they are effectively local, i.e. expressible as integrals over 4-dimensional 
Euclidean spacetime of local functions of A(x), dA(x), . . . The coefficients in these local 
functions are, of course, lattice dependent. Indeed, they must scale with the barrier height, 
A, and lattice cutoff, a -1 , to a power given by their canonical dimension. In the continuum 
limit, coefficients with negative dimensionality will tend to zero and we may therefore 
restrict attention to those with non-negative dimensionality. For the gauge variation (4.7) 
only the pseudoscalar, — needs to be considered and this functional involves only 
one relevant quantity, the dimension zero coefficient of the integral in (4.6). However, this 
coefficient generally depends in a complicated way on Aa and other dimensionless lattice 
parameters. It is expressible as an integral over the Brillouin zone and it does not agree 
with the coefficient in (4.6). But one can show that agreement is recovered in the limit, 
Aa — > 0. (This calculation was carried out in Ref.[5] where it was shown that the integral 
over the Brillouin zone develops an infrared singularity at p = in the limit A — > 0.) 
This phenomenon seems to indicate a lack of "uniformity" in the continuum limit. The 
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gauge variation of the continuum limit of T(A) does not coincide with the continuum 
limit of its gauge variation unless the secondary limit, A — > 0, is also taken. It should be 
interpreted as the lattice version of an effect that is familiar in continuum chiral theory. 
There, it is well known that the parity violating amplitude is ultraviolet convergent and 
unambiguous, as must be its gauge variation, the (consistent) chiral anomaly. On the 
other hand, this gauge variation is expressible as the difference of two formally identical, 
but ultraviolet divergent integrals, that have to be regulated carefully in order to obtain 
the correct anomaly. Ultraviolet convergent (divergent) integrals in continuum theory 
correspond to infrared singular (non-singular) integrals in lattice theory. 

In obtaining the continuum limit of the overlap amplitude we have used the approx- 
imate formula, (4.2), the leading term in an expansion of Sp in powers of p/A and ap. 
The result (4.3) is presumably valid if k±, &2, • • • , (f> are all small compared to A and a -1 . 
No condition on the magnitude of Aa is involved. However, in view of the non-uniform 
response to gauge transformations outlined above, we suspect that it would be safer to 
choose Aa to be small. 

5 Conclusions 

In this paper we have provided a set of rules for computing an overlap amplitude, order 
by order, in weak field aproximation. We have shown that this amplitude can be inter- 
preted as a lattice regularization of the vacuum amplitude for a 4- dimensional Euclidean 
continuum theory of chiral fermions coupled to background gauge and Higgs fields. 

We find that the most efficient approach is through the use of time-dependent per- 
turbation theory. The 4-dimensional Euclidean lattice is embedded in 4+1-dimensional 
Minkowski space with a continuous and unbounded time coordinate. This leads to an 
expansion of the overlap in terms of 1-loop vacuum graphs. We found it convenient in 
this work to use real time formalism but one could easily construct analogous formulae 
using imaginary time. When the background fields are slowly varying on the lattice scale 
and also on the scale of the inverse barrier height, the perturbative expressions simplify. 
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The infrared dominant term can be extracted in each order and, after integrating the time 
coordinates, put into correspondence with the continuum formula for that order. 

Although the lattice amplitude may not be itself gauge invariant, it is guaranteed that 
the continuum limit, i.e. the infrared dominant part, will be gauge invariant up to chiral 
anomalies. These would have to be compensated in the standard way in order to recover 
a fully gauge invariant theory in the continuum. 
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Appendix A: The model 

The purpose of this Appendix is to specify the details of the lattice model including the 
functions H±(p),V(p,q) used in Sec. 2 [5]. 

The dynamical variables comprise a set of Dirac fields, ipi(n,t), i = 1, . . . , N asso- 
ciated with the sites of a 4-dimensional integer lattice, n e TL A . The time coordinate is 
continuous. The Hamiltonian is bilinear and time-dependent, 

H(t) = H(n, m, t) ^(m, t) (A.l) 

n,m 

where the coefficients H(n, m, t) are 4N x 4A matrices acting in the product of Dirac and 
flavour spaces. They incorporate the couplings to external gauge and Higgs fields, A^x) 
and 4>(x), that are assumed to be smooth functions on H 4 , interpolating the lattice sites. 
The general structure is 

H(n, m, t) = H (n - m) U (n, m\A) + 5 nm M(n, t\4>) (A.2) 

where H , U and M are matrices defined as follows. 

Firstly, the gauge factor U is specified by a path ordered exponential where the path 
is chosen to be the straight line joining lattice sites n and m, 

U(n, m\A) = P (exp ji ^ dx» A°(x)T a (A.3) 

where x»(t) = tn^ + (1 - t)m^ : < t < 1. The N x N hermitian matrices T a belong to 
the algebra of the gauge group. 

Next, the mass term, M, includes the Higgs background and the barrier effect, 

75 (0 l H^ + AT c ), t>0 



M(n,t\ 



(A.4) 

[ 75 (0 i HT,-AT c ), t<0 



The N x N matrices T; incorporate Yukawa coupling parameters and define the represen- 
tation of the gauge group to which 0* belongs, i.e. 

[T i ,T a ]=i(t a ) i i Tj (A.5) 
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The "chirality" matrix, T c , is diagonal with eigenvalues ±1 corresponding to right or left 
handed flavours in the continuum limit. This matrix is required to be gauge invariant, 

[T C ,TJ=0 (A.6) 

On the other hand, since the role of Higgs fields is to connect left with right handed 
fermions, the matrices T { are required to anticommute with T c , 

{T c ,T,} = (A.7) 

Finally, to specify the hopping term H (n — ra), it is useful to employ Fourier series. 
Define the Fourier components, ip(p,t), by the lattice sum 

^(p,t) = ]T ip(n,t) e~ ipn 

n 

where pn = p^n^. These components are periodic in momenta with period 2n (in lattice 
units). The translation invariant hopping term is represented by the Fourier integral 

H (n -m) = (J^J 75(27 M C» + B(p) T c ) e^ n ~^ (A.8) 

where the integral ranges over a Brillouin zone, a cell of volume (27r) 4 in lattice units. 
The functions C M and B are real and periodic. Their detailed structure is not important 
for us, except in the infrared. We require that they have no common zeroes, apart from 
the origin. Near p = they must take the form, 

Cpip) ~ p„ + . . . , B{p) ~ r p 2 + . . . (A.9) 

where p 2 = g^p^, and r is a constant (Wilson parameter) [19]. The metric tensor, 
g^, is Euclidean and we may suppose that it is invariant with respect to one of the 
4-dimensional crystal groups. This tensor is involved also in the Dirac algebra, 

{7",7"} = 2<r (A.10) 

It is convenient to normalize the metric such that det g — 1. In choosing a metric tensor 
we are essentially choosing a crystal structure for the lattice. The particular choice is 

19 



presumably not very important as regards infrared behaviour, although we should insist 
that the invariance group of C M and B be large enough to enforce the structure (A. 9) near 
p = 0. We also assume that the crystal symmetry includes reflections with the barrier 
term transforming as a pseudoscalar. This will ensure that quantities such as $ + — $_ 
transform as pseudoscalars under space reflections of the background fields when A -1 a 
and lattice effects are ignored i.e., in the infrared regime. 
Under gauge transformations, 

A(x) -> A\x) = e i6{x) (A(x)+id) e~ i6{x) 
4>{x) ■ T -> (P\x) ■ T = e ie{x) 4>{x) ■ T e~ ie{ - x) 

the coefficient matrices (A. 2) are clearly covariant, 

H(n,m,t) -> e w{n) H(n,m,t) e~ lQ(m) 

since 6{x) = 9 a (x)T a commutes with T c . This guarantees the formula (3.2). (In the text 
we have used the collective notation, A, to represent both gauge and Higgs fields.) 

In developing the perturbative formulae of Sec. 2 we have assumed that the Higgs field 
is constant and incorporated it into the free Hamiltonian. The free Hamiltonians, Hq±, 
are therefore defined by the 1-body expressions, 

H±(p) = 75 (i Y C(p) + • T + (B(p) ± A)T C ) (A.ll) 

The vertex functions are obtained by expanding (A. 3) in powers of the gauge field. One 
obtains, 

V(p,p) = A°(k) (2tt) 4 5 2n {-p + q + k)J Q 1 dt H(p - tkr T a + 

+\ J (I 1 )' (f^) 4 A-(h) Al{k 2 ) (2vr) 4 5 2n (-P + q + k 1 + k 2 ). 

■ jf 1 dhdh H{p - tiki - hk^^h - t 2 )T a T p + 9(t 2 - ti)T p T a ) 

+ ... (A.12) 

where H{p) ,tJi = dH±{p)/dp lx , etc. The periodic delta function is defined by the lattice 
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sum, 

(27T) 4 MP) = E ^ 

n 

= J2 ( 2?r ) 4 5 *(p + 2nn ) 

n 

The momentum integrals in (A. 12) are over R 4 but, since A^x) is assumed to be slowly 
varying, its Fourier transform A^{k) is concentrated around k — 0. 

The only part of (A. 12) that is needed in the continuum limit is the first term near 
p = q = 0, 

V(p, q) = -i 757" A a (p-q)T a + ... (A.13) 
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Appendix B: Free fermions 

The purpose of this Appendix is to examine the spectrum of free fermion states and derive 
expressions for the propagators, Sp, S + and S-. 

The square of the 1-body Hamiltonian (A. 11) is proportional to the unit Dirac matrix, 

H±( P ) 2 = (T C,C U + (0 • Tf + (B± A) 2 (B.l) 

since T c 2 = 1 and {T i: T c } = 0. Choose a set of 2N orthonormal spinors x{ a ) sucn that [7] 

7sT c x(^) = X(<r), a = \,...,2N 

(<P-Tf X (a) = ml X (a) (B.2) 

and define the eigenspinors of H±(p), 

^2u±(u± + B ± A) 

I \ U ± ~ ^±(g) ^ /D 

v±(p,cr) = - === x(g) (B.3) 

/2w±(w± - 5 =f A) 



where o;±(p, <r) is given by the positive square root, 



w±(p, *) = v/CW 2 + + (5(p) ± A)2 (B.4) 

There are no zero modes, even in the absence of Higgs fields (m a = 0). The positive 
and negative energy eigenspinors, u + ,v + of H + {p) comprise a complete, orthonormal set. 
Likewise for u- and v~. The two sets are related by a unitary transformation, 

it_ = u + cos (3 — v + sin (3 

V- = u + sin (3 + v + cos (3 (B.5) 
where the angle (3(p, a) can be chosen to lie in the interval (0, 7r/2). It is given by 



/ cj + — B — A lu- — B + A u + + B + A U- + B-A 
cos/5 = W hi (B.6) 

where the roots are non-negative. 
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Near p = the energies can be expanded, 



u± = (A±B), 



1 + 



C 2 + m 2 



\ (A ± B) 2 
. , „ 1 C 2 + m 2 



2 p 2 + m 2 

A±rp + ^rr 



+ ... 



for p,m <C A. In this approximation the eigenspinors (B.3) become 

,2 \ 



M + = 



Ibiip 1 + <t> • T) p 2 + m 2 
, + 2A 8A^ + 



X 



(B.7) 



Vp 2 + 



2A 



+ ••• X 



f+ = 



l<o(ii> + <P • T) y/p 2 + m 2 \ 
+ — ^ — + ••• X 



Vp 2 + m 2 



2A 



(, Ibiif 1 + 4> ■ T) p 2 + m 
V- = 1 — Tr~m — 



and f3 approaches tt/2, 



2A 



8A 2 



+ ••• X 



COS p = : h . . . 

A 



(B-8) 



(B.9) 



To this order the quantities (B.8), (B.9) depend on the barrier height, A, but not on 
detailed lattice structure such as the Wilson parameter, r. In the next order such details 
would begin to appear. 

The long wavelength approximation to the propagators S± can be recovered from 
(B.7), (B.8) 

S±(p, t-t') = 6(t - t') J2 u ± u ± e- <w±(t - t ' ) - 

-e(t'-t) J2 v ± v ± e^-v 

a 



;b.io) 



These functions are regular at p — 0. 
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The propagator Sf is not regular at p = <p = 0. To obtain it one must consider the 3! 
possible orderings of the time coordinates, t, t' and 0, 

t(V(p,0 #z,O f ) = e(t-t')^( P ,t)^(q,ty- 

-e(t'-t)^(q,t')^(p,t) 

= e(t-t')9(t')ij + (p,t)M<i,ty 

+9(t) 9(-t') rl> + (p,t) V-(9,O t 
+0(_ f ) £( t _ t ') ^_(p >t ) ^_(g,t') f 

-6(-1?)e(l?-t)1;4q,lf)iil,-(p,t) (B.ll) 

where 

^ ± (p,t) = e _itH±(p) ^(p) (B.12) 

The matrix element of (B.ll) between free fermion ground states involves the polarization 
sums, 

<±^m = (2^MP-.)E M+(p '^- ( r )t 

(+|-> a c °s /5(p,<r) 

(+h) a cos (3{p,a) 

which are obtained by elementary considerations (see Appendix A of Ref. [7]) using the 
plane wave expansions, 

CT 

together with canonical anticommutation rules and the ground state definitions 

&±(p,g)|±) = d ± (p,<7)|±) = 

On substituting from (B.12) and (B.13) into the ground state matrix element of (B.ll), 
and using (B.5) to make the time-dependence explicit, one obtains the expression 
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S F (p,t,t') = 

= ~ *W) e ~ itW+u + («+cos/? - 4 sin/? e"*^) 

+9(t)9{-t') e~ ituJ + u+uL e tt,UJ - 

+9(-t)9(t - t') (e~ itU} - cos/? u_ + e^- sin/? u_) ul e^" 
-0(t' - i)0(i) (e- ito + sin/? m+ + e ifcJ+ cos/? u + ) uj. e- l '' w + 
-9{t')9{-t) e ltuJ - v.v\ e~* u+ 

-9(-t')9(t' - t) e ltuJ - v. (-wlsin/3 e^" + vl cos /? e^"-)] (B.14) 

This function has the expected discontinuity at i = t' but is continuous, as it should be, 
at t — and t' = 0. Its low momentum behaviour is dominated by the pole at cos/? = 0. 
It occurs at p = </> = and we can expand around this point using the formulae (B.8), 
(B.9). The result is 

Qj n t f >) - 1 + 7sTc 1 ^ a P -*A(I*M*'I) 
+ -e(t-f) e- iA l*-*'l 

+2 75 T c {(#(t - i')0(O + 9{-t')9{t' - t))e^ A l'l sin At' 

+ (fl(t' - t)6(t) + 0(-t)0(t - t'))e- im sin At} 
+ terms of order (p, </>) (B.15) 

The poles of Sf are crucial to the continuum limit since they control the infrared 
singularities. It is necessary, therefore, to establish conditions under which the pole at 
p = (ft = is unique. According to the general formula (B.14) the poles of Sf correspond 
to the zeroes of the function cos/? defined by (B.6). In this formula the square roots are 
non-negative and both of them must vanish to give a zero of cos /?. There appear to be 
two possibilities, 

lu+ = B + A, ou- = -B + A (B.16) 

and 

u+ = -B-A, tu_ = B-A (B.17) 
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The alternative (B.17) can be excluded immediately because it implies uo + + uj- = — 2 A 
which contradicts the positivity of u + and c<j_. The alternative (B.16) is possible only if 
±B + A > 0, i.e. if B 2 < A 2 . From (B.4) one sees that (B.16) implies C 2 + m 2 = 0. 
Hence, the zeroes of cos (3 occur at isolated points defined by 

C^ip) = 0, = and B(p) 2 < A 2 (B.18) 

The origin, p = 0, is certainly one such point in view of the equations (A. 9). The vector 
function, C M (p), certainly has other zeroes. This is implied by the Poincare-Hopf theorem 
since C M is defined on a torus [20]. In order that cos (3 should not vanish at these other 
points we have only to ensure that B(p) 2 > A 2 at such points. In other words, the zero 
of cos (3 at p = = is unique if B{p) and C^{p) are chosen so as to have no common 
zero, apart from the origin, and A is smaller than \B(p)\ at all the other zeroes of C^p). 
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Appendix C: Time— independent formalism 

In previous work on the overlap prescription we used time-independent perturbation 
theory and a different phase convention. The purpose of this appendix is to clarify the 
relation between the two formalisms. 

With the time-independent approach one constructs the two ground states \A+) and 
\A—) directly by solving the Schroedinger equations (2.9), 

H ± (A)\A±) = \A±)E ± (A) (C.l) 
or, rather, the equivalent integral equations 

\A±) = \±)a ± + G ± (V - E±)\A±) (C.2) 
where |±) denotes the free fermion ground states (2.5). The operators G± are defined by 

« ± = (C3) 
with the understanding that G±\±) = 0. Iteration of (C.2) leads to the formal solution 

\A±) = a ± (l - G±(V - £±)) _1 |±) (C.4) 

where the numerical factors a± = (±\A±) are determined, up to a phase, by requiring 
that the states \A±) be normalized. The energies E± are determined self-consistently 
from (C.l), 

(±\H ± (A)\A±) 
E± ~ (±\A±) 

= {±\v(l-G±(V-E±j) X \±) (C.5) 

The method sketched here is straightforward and practical, at least in the lowest 
orders. However, it is less efficient than the time-dependent method discussed in the 
text. This is mainly because it requires the calculation of the subsidiary quantities, a± 
and E±, as well as the quantity of interest, (A + \ A—). In addition, it obscures the fact 
that T(A) is expressible in terms of connected graphs. 
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A possible advantage of the time-independent method is the very simple phase con- 
vention it allows. The Brillouin-Wigner convention is expressed in the requirement that 
the numerical factors, a±, should be real and positive for all values of the external fields, 
i.e. 

(±\A±) BW > (C.6) 

This makes it easy to compute the angles &±(9,A) induced by gauge transformations. 
For example, to first order in 9, the formula (3.6) reduces to 

= Re(+\F e (l-G + (V-E + )y l \+) (C.7) 

This was used to compute the chiral anomaly in Refs.[5,ll]. 

To find the relative phase between the B-W ground state and the one used in the 
main text one must consider the definitions (2.8). These imply, in particular, 

lim (+\in+) e-* +/i£ = (+\A+) 

= (+\A+)bw e ifi+{A) (C8) 

so that f3 + can be obtained from the regular part of the connected amplitude, in the limit 
e -> 0, 

P+(A) = lim Im ( (+|m+) con - ^± ) (C.9) 

£^0 ^ IS J 

where E+ is itself given by the limit 

E + {A) = \im(ie (+\m+) con ) (CIO) 

The contribution of order N to the connected amplitude is expressed as a 1-loop integral 
constructed with the propagator S+, viz. 

(-l) N r0 



(+\in+) con = E / dh...dt N {+\T(y(t 1 )...V(t N ))\+) t 

1 V . J — OO 

= -E f d tl ...dt N e^ + - + ^. 



N 

I (fr) •••(^) tr(v(p 1 ,p 2 )S + (P2,t 1 -t 2 )...V(p N ,p 1 )S4Put N -t 1 )) 

(C.11) 
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It is straightforward to integrate the time coordinates in the terms of this series using the 
explicit formula for S+, 

S + (p,t-t>) = J2(e(t-t')u + (p,a)u + (p,ay e-^W^)- 
-9(1? - t)v+(p, <r)v+(p, a) f ef+frX*-* 7 )) 

One finds, in the limit e — > 0, 

£ + (A) = (+|V|+> + l(+|VG + V|+> + 

+^<+|^G+yG+y|+) - \(+\V\+)(+\VG 2 + V\+) + ... (C.12) 

P+(A) = ~lm (+\VG 2 + VG + V\+) + ... (C.13) 

In the same fashion one can eliminate the time integrations from the terms of the series 
(3.10) for $+ to obtain 

$+(6, A) = (+\F \+) + ±{+\(FoG+V + VG + Fe)\+) + 

+\{+\{F e G + VG + V + VG + F e G + V + VG+VG+Fg)\+) 

-l(MV\+)(M(FoGlV + VGlF e )\+)- 1 -(MFe\+)(+\VGlV\+) 

+ 1 -(+\([Fe, V]G\V - VG 2 + [F e , V]) |+) + . . . (C.14) 

To the same order the Brillouin-Wigner formula (C.7) gives 

$ + {e,A) BW = (+\F e \+) + ±{+\(F e G + V + VG + F e )\+) + 
+ l -{+\(F e G + VG + V + VG + VG + F e )\+) - 
-\(+\V\+){+\(FeGlV + VGlF e )\+) + ... (C.15) 

These phases are related by 

$+(9, A) = A) BW - 5 e P + (A) (C.16) 

where f3 + is given by (C.13). This can be verified using the formula (3.2) or 

5 V = i[F g ,H o+ + V] 
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According to (C.16) the chiral anomalies are related by 

and we must consider what happens to the functional, f3 + — (3_, in the continuum limit. 
Like $±, the angles (3± do not involve any infrared singularities. In the continuum ap- 
proximation they must be local, i.e. expressible as integrals over 4-dimensional spacetime 
of local functions of the gauge field and its derivatives, with lattice dependent coeffi- 
cients. More particularly, the difference, (3 + — must involve a pseudoscalar density. 
For example, 

(3 + -(3-~^J d 4 x g o° e^ v d K A x d p A, d a A v + . . . 

There is no candidate of dimension 4. This means that in the continuum limit, k/A — > 0, 
this functional becomes vanishingly small: the continuum theory chiral anomalies are 
unaffected by the phase conventions. 
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